The Impact of Inhomogeneous Reionization on 
Cosmic Microwave Background Anisotropy 



00 



Lloyd Knox^, Roman Scoccimarro^, and Scott Dodelson^ 
^ Canadian Institute for Theoretical Astrophysics, Toronto, ON MBS SH8, CANADA 
^ NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 
(February 1, 2008) 

It is likely that the reionization of the Universe did not occur homogeneously. Using a model 
that associates ionized patches with overdense regions, we find that the resulting cosmic microwave 
background (CMB) anisotropy power spectrum peaks at angular scales corresponding to the extent 
of the ionized regions, and has a width that reflects the correlations between them. There is consid- 
erable uncertainty in the amplitude. Neglecting inhomogeneous reionization in the determination 
of cosmological parameters from high resolution CMB maps may cause significant systematic error. 
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Introduction. Observations of cosmic microwave back- 
ground (CMB) anisotropy are providing strong con- 
straints on theories of cosmological structure formation. 
Planned CMB observations can potentially provide con- 
straints on the parameters of these theories at the percent 
level (|,|. One of the reasons the CMB is such a wonder- 
ful probe of cosmological parameters is that predictions 
for a given theory can be made with great precision 
This is because linear perturbation theory is an excellent 
approximation over the relevant length scales at z ~ 1100 
when the CMB was last tightly coupled to matter. 

However, new anisotropics can be generated at late 
times by the non-linear process of reionization as the 
CMB photons are brought back into contact with matter 
via Thomson scattering. We know this reionization took 
place before redshift z ~ 5 because spectra of distant 
quasars do not show a continuum of absorption by neu- 
tral hydrogen. It is unlikely that reionization occurred at 
z > 60 for in that case the level of anisotropies on degree 
scales would be significantly lower than on ten degree 
scales, when actually just the opposite is true Q. 

Reionization affects the CMB in several ways. First, 
the fluctuation power on small angular scales is damped 
by a factor e"^"^ where r is the optical depth back to last 
scattering. Second, as Sunyaev and Kaiser j|] pointed 
out, secondary anisotropies are generated on large scales 
due to the Doppler effect when photons scatter off moving 
free electrons. They also noted that the effect is strongly 
suppressed on small scales because photons get nearly op- 
posite Doppler shifts on different sides of a density peak, 
a consequence of potential flows generated by gravita- 
tional instability. Both of these effects — the damping and 
the Sunyaev-Kaiser (SK) effect — are linear and therefore 
included in standard Boltzmann codes 

Here we focus on the effects of inhomogeneities in 
the ionization field We show that inhomogeneous 

reionization (IHR) restores the SK effect at small scales 
due to modulation of the velocity field with the spatial 
variation of the ionization fraction, an effect analogous to 



the modulation by the density field in the homogeneous 
case at second order in perturbation theory |^ . Since the 
typical patch size is much smaller than the scale of varia- 
tion of the velocity field, the Doppler effect contributions 
come only from regions where the flow is coherent. The 
resulting anisotropy pattern depends sensitively on the 
model of reionization. Thus, the signature of IHR in the 
CMB provides an indirect window on a poorly probed 
but interesting era in which the first stars and galaxies 
formed. On the other hand, if these anisotropies are large 
enough, they can create systematic errors in the deter- 
mination of cosmological parameters. So, after calculat- 
ing the anisotropies, we compute the magnitude of these 
systematic errors, i.e. how much a given cosmological 
parameter might be misestimated if IHR is ignored. 
Anisotropy Spectrum: Formal Solution. The perturba- 
tion to the photon temperature distribution function, 
A = 5T/T, is governed by the Boltzmann equation, 
which at the late times of interest is 



5A(x, 77, ^)/dri = Heara-j ■ v(x, 77) 



(1) 



where a = (1 + z)^^ is the scale factor of the universe 
normalized to unity today, rj = J dt/ a is the conformal 
time, Ue is the free electron density, ctt is the Thomson 
cross-section, 7 is the direction of the photon momentum, 
and V is the electron velocity. Henceforth we work in 
units where the conformal time today is unity, rjo ~ I; 
and assume a flat, matter dominated universe, a = rj^. 
The solution to Eq. |l| for the photon perturbations here 
and now (at x = xq and ry = 1; Ao(7) = A(xo, ?7o, 7) ) is 



Ao(7) = 4.06 X 10"^ flbh 



^ drj / N 

^ Xe x,?7 7-t;(x , 2 



where, v(x, 77) = —2r] t>(x) (Mpc/h), Ue — 
•he Xe(x, 77)0"^, and neCTT^o = 0.122ilbh, with fie the 
present mean electron density, fit the baryon density in 
units of the critical density and h the Hubble constant in 
units of 100 km s~^ Mpc~^. Both Xe and v are evaluated 
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at X = Xo — 7(1 — 'y); the position at time rj oi a. photon 
with direction vector 7 incident on us today. Note that 
the integral includes all contributions starting from rji, a 
time far after standard recombination at 77 ~ 0.03, but 
before reionization occurs. 

We focus on predicting the two-point correlation func- 
tion, C{9) = (Ao(7)Ao(7')) where cos 9 = 7-7' and the 
angular brackets indicate an average over all locations 
xq . Using Eq. |2| we obtain the correlation function due 
to reionization: 

It is therefore determined by an integral over adjacent 
lines of sight of the four-point correlation function 



C{r],ri',9) = 7i7j (xex'^ViVj) 



(4) 



where, for example, Xg = Xe (x',77'). Different models 
of reionization will lead to different four-point functions. 
We now explore several possibilities. 
Homogeneous Reionization. The simplest possibility is 
that reionization takes places homogeneously. While 
not particularly plausible, it is a useful limiting case 
for demonstrating the SK cancellation. In homogeneous 
reionization (HR), Xe = Xeif]), where Xe is fraction of the 
Universe that is ionized. Thus we have: 



(5) 



and we take Xe to be zero at times earlier than Zi -\- Szi, 
unity at times later than Zi and to increase linearly with 
redshift during the transition. The velocity correlation 
function can be written as 



a(77,77',0) = *xcos0+(*||-*_ 



(r • 7) (r • 7') 



, (6) 



where r = x— x'. The correlation functions of the velocity 
components parallel and perpendicular to r are 



'f|iW = 



P{k) 



Jiikr) 



kr 



m Jijkr) ^3^^ 
k'^ kr ' 



(7) 
(8) 



where P{k) is the matter power spectrum, throughout 
assumed to be that of the standard cold dark matter 
theory [|o|. Fig. 1 shows Ci,(ry,7/,0) as a function of 77 
for fixed 77' = 0.188. is strongly peaked at equal times, 
where the two lines of sight share the same velocity. On 
the other hand, at large 177 — 77'!, Cy becomes negative 
due to infall from opposite sides into an overdense region. 
The integration, in Eq. ^ of this oscillatory function with 
the slowly- varying x^x'^, leads to a cancellation of Doppler 
effects, as described by Kaiser Q in Fourier space. 
Uncorrelated Inhomogeneous Reionization. The SK can- 
cellation can be avoided if the velocity two-point func- 
tion Cy is modulated by the spatial dependence of Xe- 



To demonstrate this effect, we adopt a toy model of the 
reionization process: independent sources turn on ran- 
domly and instantaneously ionize a sphere with comov- 
ing radius R and volume Vr, which then remain ionized. 
The resulting four-point function can be written as: 



C-(77,7;',0) =^(77,77',^) [C::ifj,rj',9)+x,x', 



(9) 



Since x (where x = ^~Xe) is equal to the probability that 
no ionizing source is within volume Vr, the probability 
that no source is within R of either x or x' is (at equal 
times) ixx') = = XXX'^^'^'^^^"' where for 

r < 2R, Vir) = Vr [I - (3/4) (r/i?) + (1/16) (r/i?)^] 
is the overlap volume of the spheres centered on x and 
x'. Therefore, with the proper generalization to unequal 
times, we have: 

^ (x^X) - ^eK = ixx') - xx' 
= xx' 



1 

2 0.1 

°. 0.01 

C 10-3 

°. 10-= 
— 10-' 
10-" 



x(min(7/,77')) 


-V(r)/VR 


- 1 




; y.-^ — y 












" , , , , 1 , 


1 , , 




18 0,185 


0.19 


0.1 



(10) 

(11) 



FIG. 1. The mean ionization fraction Xe for HR {zi = 26, 
Szi = 4), the velocity two-point function C„ (77, 0.188, 0) for 
the SCDM (standard cold dark matter) model with as = 1.2 
where as is the rms of mass fluctuations in spheres of radius 
Sh~^ Mpc today, and the ionization fraction two-point func- 
tion for uncorrelated C"^ and correlated bubbles CJ° . Dashed 
lines denote negative values. 

The resulting two-point function, C^'^, is sharply 
peaked at small r w [77 — ri'\ (see Fig. 1, where R — 1 
Mpc), therefore in Eq. |, ^"{9) cx Cy{7j, 77, 9). Thus, the 
SK cancellation of homogeneous reionization is avoided; 
IHR modulates the velocity field, accessing only the re- 
gion where the velocity field is highly coherent. Note 
though that the modulation is not particularly effective 
since the ionization radius R is typically very small, and 
thus only a tiny region of rj contributes. Therefore, as 
Gruzinov and Hu (GH) Q have pointed out, a model 
of this type necessarily produces anisotropics which have 
amplitude proportional to R. 

Figure 2 shows the anisotropy spectrum arising from 
the two-point function of Eq. ^ assuming Zi = 26 and 
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5zi = 4. The GH results for the same parameter choices 
would each have peaks 30% lower in / and two times 
larger in amplitude due to a slower fall off of C^'^ with r 
and one less factor of {1 — Xe) in their version of Eq. 
Despite these differences, they correctly identified the 
qualitative features of this type of IHR: (i) amplitude 
proportional to R; (ii) white noise (C™^ oc constant) on 
large scales due to lack of correlations at r ^ i?; and (iii) 
peak at I ^ {1 — r]r)/R where reionization takes place 
at rir- We now show that the first two of these features 
break down when more realistic models are considered. 




103 10" 
multipole moment, I 

FIG. 2. The angular power spectrum, Ci, defined by 
C{e) = Y.A'^^ + l)Ci/(47r)P;(cos6'). The thick solid line on 
top is the linear spectrum. The rest are different approxi- 
mations to the IHR anisotropies, all assuming 57;, = 0.1 and 
h = 0.5. The curves peaked at I = 4000 and I = 20000 are 
from the uncorrelated _R = 5 Mpc and R = 1 Mpc models 
respectively. The oscillations are ringing due to the sharp 
cutoff of V{r). The dashed curves show the model of Eq. [13 
for Zi — 26 (bottom) and Zi — 31 (top). The remainig solid 
lines include velocity-ionization field cross-correlations. The 
correlated model curves do damp beyond £ ~ 1/R > 30000. 

Correlated Inhomogeneous Reionization. Because the 
ionizing sources in the model of Eq. |ll| and that devised 
by GH are independent of each other, correlations in the 
ionization field only exist on scales comparable to, or 
smaller than, the ionization radius R. However, over- 
dense regions, such as galaxies and clusters of galaxies, 
are observed to be highly correlated with each other. If 
overdense regions are the sources of the ionizing radiation 
(which they are very likely to be) then their clustering 
will lead to long-range correlations in the ionization field 
and therefore affect C™^(6'). 

The physical process we imagine is that the mass in any 
region where the linear theory density contrast smoothed 
on a comoving scale i?c, Sr^, exceeds a critical threshold 
(5c — 1.69), collapses, forms stars or quasars which then 
ionize a region with size R = E^I^R^. The efficiency, 
i?, is the ratio of the volume of the ionized regions to 
collapsed regions. From this scenario we expect ]To| , p^ 



= (xcx'c) = OLol / dS / d8' Pr,{S,S') 



aa 



27r CTfl,^ (77)0-^^(77') 



-exp 



(12) 



, (13) 



where erfc(x) denotes the complementary error function, 
a = 2E{1 — Xe) (see below), Pr^{S) and Pr^{6,6') are the 
one and two-point distributions of the smoothed field 6r^ , 
respectively; the smoothed correlation function is given 
by (56') = Sfl^(r), aR^ = Sfl;^(0), and 1/(77) = Sc/ctrXii)- 

Despite a number of studies of the details of reion- 
ization [pl]~p^, we are still far from a satisfactory un- 
derstanding. Therefore, we take two simple models for 
E, which can vary by many orders of magnitude (see 
e.g. table 2 in ||l^). The requirement that objects 
more massive than Mc{z) ~ 1O^M0[1O/(1 + z)fl'^ are 
necessary for reionization leads to a smoothing scale of 
Rc = 0.22r/ Mpc p^. Our two choices for the efficiency 
are £; = 7 X 10^ vf" resulting in z, = 26 and E = 114 
(the "middle of the road" model in (l^), which gives 
Zi = 31. Note that current scenarios for reionization 
generally lead to i? ^ 1, i.e. the sources tend to ion- 
ize a region much larger than their typical volume (see 
however |^ for a different view). 

Fig. 1 shows the results of Eq. |l^ for the Zi = 26 model. 
The clustering of overdense regions naturally leads to 
a much wider C^°, which falls off only as a power-law, 
thus more efficiently modulating C^. For this reason, 
the correlated models in Fig. 2 (dashed lines) show a 
much wider distribution of power, and the white noise 
regime is only reached at scales much larger than the 
patch size. Over the range of I plotted, the power spec- 
trum of the uncorrelated models with the same values 
of R as in the peaks model would be orders of magni- 
tude smaller. The dashed lines in Fig. 2 are the result 
of assuming (xex'^vv') = (xex'^) (vv') and then applying 
Eq. 13 for {xex'c). For the solid lines, we use an expres- 
sion for (xex'^vv') similar to that for {Xf,x'^) in Eq. 13 
but with an integral over the joint probability distribu- 
tion of (5, (5', w, v' This leads to a modification of the 
white noise behavior at large scales. Cross-correlation 
between HR and IHR leads to higher-order terms such 
as Xe{x'^Svv'). These are of comparable magnitude to 
{xex'^vv') and will be discussed in |l^. Note that the 
fiatness of the power spectra in Fig. 2 above £ ~ 3000 re- 
fiects the near scale-invariance of the matter power spec- 
trum over the relevant length scales. 

Impact on Parameter Determination. The anisotropies 
produced during IHR may significantly affect the 
exquisite parameter determination anticipated from fu- 
ture CMB observations. Suppose a multi-dimensional fit 
is performed on the data to extract a set of parameters 
{pi\. If the fit assumes Ci{{pi]) that ignores the con- 
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tribution from C™^, then each parameter will be incor- 
rectly estimated by an amount 



dpj 



qIHR 



(14) 



where the weights wi are the inverse of the squares of 
the errors expected on C/'s (the errors are the sum of 
the errors due to sample variance and those due to noise 
§,1); is the Fisher matrix J2i wi{dCi/dp^){dCi/dp,). 
Table 1 shows the ratio of this systematic offset to the 
statistical uncertainty for Planck . 
Sources of Uncertainty. Even within the context of our 
peaks model, Eq. |l^ is not exact. Since the ionized re- 
gion is larger than the region that collapsed (i.e., E > 1) 
a more rigorous approach would be to calculate the prob- 
ability of two points both being within a distance R of 
one or more peaks. Instead, we have calculated the prob- 
ability that two points are both within a peak and then 
scaled the resulting two-point function by EE' . We ex- 
pect this approximation to be best for r > R, which 
covers the range of relevance for the curves in Fig. 2. Be- 
cause we have not taken the more rigorous approach, we 
must include suppression factors to correct for overcount- 
ing overlapping regions. These are the factors of (1 — Xe), 
motivated by their appearance in our fully tractable un- 
correlated model. Removing them boosts the power spec- 
trum by a factor of about 5. 

Uncertainty in E has a much milder effect on uncer- 
tainty in C™^ than one might expect from a quick glance 



at Eq. 13 



The reason is that increasing E causes reion- 
ization to occur earlier and therefore decreases the expo- 
nential factor in Eq. |l^. We find that decreasing E hy a 
factor of 10 decreases C™^ by a factor of two. 

Our peaks model, and the calculations which lead to es- 
timates of its parameters, are highly idealized. Although 
more sophisticated studies will eventually alter the de- 
tails of our results, we believe significant power at scales 
larger than the patch size is a natural consequence of 
cosmic structure formation via gravitational instability. 
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Zi = 31 
Zi = 26 


0.75 
0.34 


0.65 
0.29 


0.01 
0.006 


2.25 
1.01 


0.20 
0.09 


-1.02 
-0.46 


0.94 
0.42 



TABLE I. Ratio of systematic offset to statistical uncer- 
tainty expected for Planck |l^ from the two "peaks" models 
in Fig. 2. Since IHR affects very small scales, MAP |l|] 
is much less sensitive to this offset. Variables are, from left 
to right, the primordial power spectrum index, the ratio of 
tensor to scalar contributions to C2, the optical depth and 
the contribution to from (b)aryons, (curv)ature, cold dark 
(m)atter and massive neutrinos. The power spectrum deriva- 
tives were evaluated at the parameter values of the SCDM 
model in [^. This does not include polarization information; 
however, since polarization is sourced by the quadrupole mo- 
ment, not velocity, the polarization IHR power spectrum will 
be down by roug hly a factor of (10"VlO"^)^ from the IHR 
temperature power spectrum. 
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